New Quantum Bounds for Inequalities involving Marginal Expectations 
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We review, correct, and develop an algorithm to determine arbitrary Quantum Bounds based on 
the seminal work of Tsirelson [Lett. Math. Phys. 4, 93 (1980)]. The vast potential of this algorithm 
is freshly demonstrated by deriving both new number-valued Quantum Bounds, as well as identifying 
a new class of function-valued Quantum Bounds. Those results facilitate an 8-dimensional Volume 
Analysis of Quantum Mechanics which extends the work of Cabello [PRA 72 (2005)]. Finally we 
contrast these bounds to those defined by the first-level criteria of Navascues et al [NJP 10 (2008)], 
proving our function- valued quantum bounds to be relatively more complete. 
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I. INTRODUCTION 

Quantum Mechanics (QM) is an inextricably nonlo 
cal |l| theory: it is fundamentally incompatible with Lo 
cal Hidden Variable Models (LHVM). This is perha 
best evidenced by the violation of Bell inequalities 
3] , although other manifestations of nonlocality are also 
known [4J47|. On the other hand it is now appreci- 
ated that QM is less than maximally nonlocal, in the 
sense that causality-respecting, i.e. no-signaling-fastcr- 
than-light (NOSIG) frameworks allow for violation of the 
Bell inequalities even beyond the Quantum limit 
Many recent works have offered information-theoretic ex- 
planations to account for the discrepancy between those 
statistics a-priori permitted by General Relativity and 
those a-posteriori limited by quantum theory 
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In this work we investigate the distinctions between 
LHVM, QM, and NOSIG through gedankenexperiments 
consisting of parties (A, B,C . . .), each having access to 
measurement apparatuses (Aq, A\ . . . Bq, B\ . . .) capable 
of finite measurement outcomes. Here we consider only 
binary measurement outcomes, which we take to be ±1. 
Practically this amounts to spin measurements on physi- 
cally separated qubits, where the Aq is a spin measure rel- 
ative to angle 0a o and so forth. The statistical quantities 
of interest are the marginal probabilities P(Aq = —1), 
P(Aq = +1) for all the measurement apparatuses for all 
the parties, as well as the joint coincidence counts such 
as P{Ai = Bj) which we can equivalently express more 
generally as P ((Ai x Bj x Ck x ...) = 1) for all i, j, k, ... 
Here we follow the conventional notation of expectation 
values such that (A4) = P(A l = 1) - P(Ai = -1) 
and (Ai = Bj) = (Ai ■ Bj). We'll focus on inequal- 
ities (bounds) pertaining to linear combinations of the 
expectation values, such as 

(A ■ B ) + (A ■ B 1 ) + (A 1 ■ B ) - {A 1 ■ B x ) < 7 (1) 

LHVM posits that all measurements merely reveal pre- 
established values. One can demonstrate (1(3, E3] that 
all statistics consistent with LHVM "fit" within Eq. (|TJ) 



upper-bounded by 7 = 2. One can also prove 0, [Io| 
that NOSIG, in contrast, permits Eq. ([1]) to go all the 
way up to 7 = 4. The quantum bound for Eq. ([TJ has 
been proven [l5| to be exactly 7 = 2\/2, in between the 
bounds of LHVM and NOSIG. We'll refer to the minimal 
binary and dichotomic scenario implicitly referenced in 
Eq. §Q as the EPR-Bcll scenario [IH1] ; it is the most 
fundamental setup which allows for nonlocality, defined 
as an inequivalence between the statistics of LHVM and 
NOSIG. 

To investigate the distinctions between LHVM, QM, 
and NOSIG models we contrast the different limitations 
they imposed on possible gcdankcncxpcrimental statis- 
tics. Finite and complete sets of inequalities have been 
determined for both NOSIG and LHVM in many mul- 
tipartite and multichotomic scenarios IB - HI, however 
identification of the restrictions due to QM remains an 
open problem, even when limited to only the EPR-Bell 
scenario. A necessary and sufficient characterization of 
QM has been identified for the 4~Spa.ce of correlations 
within the EPR-Bell scenario, in which all marginal ex- 
pectation values, e.g. (^4^), are taken to be identically 
zero. We refer to those conditions as the TLM*- 4 - 1 in- 
equalities, after Tsirelson j^, Landau [Hj], and Masanes 
[22J) who each derived an equivalent form of them in- 
dependently. Efforts to address the full 8-Space of both 
marginal and joint probabilities have been less developed. 



Tsirelson shared a theorem on the matter in 1980 [15 1 
followed by the hierarchy of semi-definitc-programming 
tests of Navascues, Pironio, and Acm (NPA) in 2007 
Hill, which converge in their infinite limit to a com- 



plete characterization of QM. Application of the NPA 
algorithm to the first-level of the hierarchy yielded the 
NPA 1 - 8 ' inequalities, the first (and heretofore only) known 
inequalities in 8-Space. The NPA (8) inequalities exactly 
identify statistical behaviors consistent with the princi- 
ple of Macroscopic Locality [25| . Those conditions have 
been used since as a 'gold standard' to approximately 
delineate QM (8) , such as in Ref. J26|, but one should 
not mistake them for complete. Here we have introduced 
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our superscript notation to indicate the probability space 
acted on by a given set of inequalities [27 1 . 

In this paper we obtain a new characterization of 
QM^ 8 - 1 distinct from, and more complete than, NPA^ 8 -* . 
To do so we requisitioned Theorem 2 of Ref. |15| , which 
was published without proof [28j|. In the process we 
identified a critical mathematical error which had been 
heretofore unnoticed. In this paper we construct the al- 
gorithm from first principles anew, thus providing both 
proof and correction to Tsirelson's theorem. The specific 
repaired expressions appear in Eq. ([5]) here. For peda- 
gogical purposes, as well as to facilitate a comparison 
to NPA (8) , we develop the algorithm primarily within 
the EPR-Bcll bipartite dichotomic scenario. We stress, 
however, that this algorithm naturally generalizes to any 
binary setup and is always simulated by shared qubits. 



II. THE ALGORITHM 

Consider the following highly-general question: Given 
an arbitrary 8- Space measure (Z), weighted by real pa- 
rameters Ci such that 

Z= Cl A Q + c 2 A 1 + c 3 B a + C4B1+ (2) 
c 5 A ■ B a + c 6 A 1 -B a + c 7 A ■ B 1 + c 8 A x ■ B 1 

what is the upper limit on (Z); what is (Z) Max =? 

For instance, choosing 1 = C5 = cq = C7 = — eg with 
all the other a = yields Eq. (JTJ). Determining upper 
bounds for such linear combinations of expectation val- 
ues is a relatively straightforward problem in both LHVM 
and NOSIG models, as the inequalities which perfectly 
define them in 8-Space are known [13, [21 ■ How- 

ever, since no such perfectly-defining inequalities exist 
yet for QM, a unique algorithm is called for, which we 
present here. 

First we analogize the scenario's measurement appa- 
ratuses into quantum measurement operators. The spa- 
tial separation requirement is enforced in the condition 
that the operators of the different parties Ai, Bi,... com- 
mute with each other. We choose to put the operators in 
separate Hilbert spaces [3(J 31 1. The ±1 possible mea- 



surements are embedded in the operators' eigenvalues, 
e.g. spin projectors to the x and y spin axes. Tsirelson 
proved that a 2-dimcnsional Hilbert space for each party 
was sufficient to emulate any general binary correlation 
behavior . Without loss of generality we can consider 
only the difference in angle between a party's projective 
measurements, such that we can impose a reflection sym- 
metry for dichotomic scenarios pjij ]. i.e. 



A Q = (cos(0a)<Jx + sin(6A)<J y ) ® 1 (g> 1 . 
Ai = (cos(9a)& x — sin(9i)a 9 ) <g> 1 <g> 1 . 
B a = 1 (g) (cos(8b)cTx + sm(9 B )a y ) ® 1 . 
Bi = 1 (8 (cos(8b)(7x - sm(9 B )a y ) ® 1 . 



(3) 



The quantum bound is, equivalently, the largest possi- 
ble measurement value of Z, the supremum of the largest 
eigenvalue, the global maximum of the right-most root of 
the secular equation. For the 8-Space Z defined in Eq. ((JJ 
the secular equation has four roots, each parametrized by 
only the two variables, Ba and 9b- Writing the charac- 
teristic polynomial of Z in terms of the variable m we 
identify the quantum bound for Eq. ([2]) as 



(Z) 



Ma 



sup(m) 



such that 



where 3 |u| 
/x 2 m 2 + [i 3 m - 



= H = i 

fj,4 = 



(4) 



Here we are using a corrected form of Tsirelson's compact 
notation [15| such that 



/i 2 = -(|e| 2 + |/| 2 )-2(|. 9 | 
fj.3 =4Rc[fhg + ehg] , 

(Iffl 2 



/'4 



\e\VV 



N 2 ), (5) 
2Re[feh 2 + eff] , 



e = UVC5 + vucq + UVC7 + uvc$ , 

/ = uvc$ + uvce + uvc-j + vuc$ , (6) 

g = uc\ + UC2 , h = VC3 + VC4. , 

where u & v are complex unitary variables, such that 



u = cos(9a) + i sin(^) . 
v = cos(9b) + isin(^s) ■ 



(7) 



and the over-bar notation in Eq. ([S]) and Eq. (JB]) is under- 
stood to indicate complex conjugation. 

Eq. ((4|) requires the analytically-intractable steps of 
identifying the largest solution to a quartic equation as 
well as maximizing over two variables simultaneously 
[33| . We bypass both obstacles through the use of for all 
statements, shifting the variation from u and v to only m. 
Observe that the characteristic polynomial has a definite 
positive leading coefficient (namely 1), and thus our op- 
erator's largest eigenvalue corresponds to that root of the 
polynomial at which all the derivatives are non-negative. 
As such, the quantum bound is given by the alternative 
formulation 



(Z)mox = inf(m) such that 

™4 1 ,, ™2 



V \u\ = \v\ = 1 



and 
and 
and 



Ii2m' 
4m 3 



- ^i 3 m + ^4 > 

2/12TO + [A3 > 

6to 2 + fj, 2 > 

m > 



(8) 



which is now suitable for analytic analysis. Note that 
the number of positivity constraints scales like two to the 
number of parties, and does not vary with multichotomic 
generalization. 
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TABLE I. bounds (Number- Valued and Function- Valued): TB' 4 ' is the well-known Tsirelson's bound, and QB/ 4 ' is its function- 
valued generalization. QB 3 ( - 8 ' ) surprisingly coincides with LHVMMax < NOSIGMax for all 1 < |x-| < 2. QB 3 ' 8 ^ is the most 
powerful of the bounds, alone capable of excluding by almost twice the volume excluded by all other known bounds combined, 
including NPA (8) . 
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\x\+2 V \x\ > 2 
V2x 2 + 8 V \x\ < 2 



2 V \x\ > 2 

V |x| < 2 



QB 3 (8) 



-1 



3|x| - 2 for |z| > 2 
b| +2 for b| < 2 



3|x| - 2 
b|+2 

x 2 



3x :- i0x '? a v |x| < i 



V b| > 2 

V 1 < k| < 2 j3|x|-2 Vb|>2 
4 V bl < 2 



III. NEW QUANTUM BOUNDS 

Our ambition was to study those sets of Ci where the 
LHVM and NOSIG bounds do not coincide. Wc decided 
to examine all possible unitary or zero values for the eight 
Ci = — 1, 0, 1. A naive approach would be the computa- 
tion and collection of all 3 8 = 6561 sets, but through var- 
ious physical symmetries (see Appendix B) we reduced 
those to 98 inequivalent candidates, of which only three 
yield a gap between LHVM and NOSIG, listed in Ta- 
blcUl TB*- 4 -* there corresponds to Eq. ([T]), but other two 
number-valued bounds are new to the 8-Space analysis. 

The algorithm's most powerful feature is not its yield- 
ing of number-valued quantum bounds; other numer- 
ical algorithms can do so more efficiently [HI]. The 
more significant reward of this analytical algorithm is 
that, when paired with the symbolic algebra prowess of 
Mathcmatica™, our algorithm unleashes a new class of 
function-valued quantum bounds, such as QB 1 ^ 4 \QB 2 ( - 8 \ 
and QB 3 (8) in TableU 

Interestingly QB 3 ' 8 - ) demonstrates zero non- locality 
beyond the LHVM models for \x\ > 1 despite the over- 
head permitted by NOSIGMax for Ixl < 2. This would 
allow for an "inverse Hardy" -type 35, 3(| all-or-nothing 
test of Quantum Mechanics, in that the presence of any 
non-locality would not vindicate Quantum Mechanics, 
but rather contradict it for the parameter region in ques- 
tion. 



IV. VOLUME ANALYSIS 

We can now contrast the function-valued quantum 
bounds in TableU with the known quantum bounds 
TLM (4) and NPA (8) . Both TLM (4) and NPA (8) have the 



same form, namely that for all i and j, 

|f(o,o) + f(o,i) + f(i.o) + f(i,i) - 2 f(j,fc) | < 7T 

where for TLM 1 - 4 -', as derived by Masanes [22| 
f(j,k) = arcsin(AjBfc) 

and for NPA (8) the authors found 



f(j,k) 



arcsm 



(AjB k ) - (Aj) (Bj 




(9) 



(10) 



(11) 



Recall that TLM^ 4 - 1 is complete in 4-Spacc, whereas 
NPA^ 8 ) is incomplete in 8-Space. Note that by setting 
the marginal expectations to zero in NPA (8) one recovers 
TLM (4) , and therefore TLM (4) is wholly included within 
NPA (8) . 

We quantify the tightness of a set of criteria using 
a technique called Volume Analysis, first introduced in 
Ref. Q; see therein for physical interpretations of the 
volume. As an informal estimate of the completeness 
of function- valued bounds in general we contrast the vol- 
ume of QB^ 4 ' to both the weaker TB^ 4 - 1 and the stronger, 
complete, TLM (4) in 4-Space. The 4-Spacc Volume Anal- 
sis for TB (4) and TLM (4) has already been done in Ref. 
, however the analytic value of the volume for TB*- 4 -* 
is original. We find that the function- valued QB^ 4 - 1 re- 
sults in a significant tightening compared to the number- 
valued TB (4) , in that QB^^ more than halves the excess 
volume of TB^ 4 - 1 relative to TLM (4) . 

To analyze a function- valued bound such as QB^ 4 -* we 
compose a Boolean function of x that takes on the value 
1 iff the inequality is violated, and returns otherwise. 
A point, i.e. a set of hypothetical statistics, is within 
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TABLE II. Results of 4-Space Volume Analysis. The entirety 
of 4-Space has volume 2 4 = 16, and the tabulated values are 
understood as multiples of 2 4 . Note that the all points in 
4-Space fit within NOSIG (8) . 



TB (4) 


QB/ 4 > 


TLM (4) 


LHVM (4) 


48v|-65 „ g 61 


« 0.938 


^ « 0.925 


| « 0.667 



the volume of this pseudo-quantum polytope is somewhat 
larger than w 1036 X This lower bound for QM^ is so 

9 8 

far below the upper bound of 1083.8 x ^ such that it does 
not offer much insight. Polytope construction with addi- 
tional extremal points would undoubtedly yield a better 
lower bound, however such higher-precision efforts were 
effectively unfeasible using our non-professional compu- 
tational methods. 



the bound iff the integral over x from — oo to oo is pre- 
cisely 0. We used a Monte Carlo method to estimate the 
volume, sampling many random points to determine the 
approximate fraction contained within the criteria. 

To construct an approximation for QM^ 8 -* we combined 
all available quantum bounds, namely NPA'- 8 ' 1 , QB 2 ^ 8 - ) , 
and QBg' 8 - 1 and applied them to random nonlocal points, 
i.e. statistics rejected by LHVM*- 4 -' but permitted by 
NOSIG (8) . We do not include QB^ 4 ) because its restric- 
tions are wholly contained within TLM (4) which in turn 
is contained entirely within Using this model we 

found an upper bound for the true quantum volume in 
8-Space of 1083.8 x We also originally and analyti- 
cally determined the volumes of LHVM (4) and NOSIG (8) 
within the 2 8 = 256 total volume of 8-Space. 

We found QB3' 8 - 1 to be the most powerful bound, in 

9 s 

that it alone has a volume of ss 1084 x jg, apparently 
singularly excluding 99% of the nonlocal points excluded 
by the best composite model. Further composite models 
without QB 3 ^ all had volumes larger than 1085 x ^ . We 
also found the volume of NPA (8) to be « f 085.8 x sig- 
nificantly less than the volume of either QB 2 ( - 8 ' ) or QBg*- 8 ) 
alone. We did identify some points that were forbidden 
by NPA< 8 > but not by QB 2 ( - 8 ^ nor QBg*- 8 ), however these 
exceptional statistics comprise a relatively miniscule vol- 

9 s 

ume, numerically capped by at most 0.1 x . To summa- 
rize: QBg 1 - 8 ' is significantly more restrictive than NPA (8) , 
and NPA*- 8 ' is almost entirely contained within QB 3 ^ 8 - 1 . 

It would be gratifying to have some lower bound for 
the true volume of QM*- 8 '. Towards this goal we con- 
structed a pseudo-quantum polytope defined by the 16 
extremal points of LHVM' 4 - 1 together with the quantum 
analog of the 8 nonlocal extremal points of NOSIG' 8 - 1 . By 
a variation of the convex hull method we determined that 



TABLE III. Results of 8-Space Volume Analysis. The entire 
probability space, within which No-Signaling theories are con- 
tained, has total volume 2 8 = 256, and the tabulated values 
are given as multiples of 2 8 /8!. The volume associated with 
QM 1 - 8 ' below is necessarily an overestimate, as the compos- 
ite model used for QM' 8 ^ is only a collection of incomplete 
criteria. 



NOSIG (8) 


N PA( 8 ) 


« QM (8) 


LHVM (4) 


= 1088 


w 1086 


< 1084 


= 1024 



V. CONCLUSION 

Thus we have rederived, repaired, and repurposcd 
Tsirelson's theorem (l5j into an algorithm for obtaining 
analytic quantum bounds. In our formulation it is ex- 
plicitly suitable even for the study of multipartite nonlo- 
cality, and in principle extendable also to multichotomic 
scenarios, although we restricted our focus in this pa- 
per to the 8-Space of the bipartite dichotomic EPR-Bell 
scenario. We identified new quantum bounds relevant 
to this highly-studied scenario, and then used those new 
bounds to perform an approximate volume analysis for 
QM (8) . 

This algorithm is only one path towards approaching a 
genuinely complete characterization, the convergent hier- 
archy of Navascues et al [24[ is another. Still missing, but 
eagerly anticipated, is the discovery of complete inequal- 
ities defining QM^ 8 * 1 , akin to the discovery of TLM^ 4 * 1 for 
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APPENDICES 



APPENDIX A: The Algorithm as Computer Code 



Here we present the algorithm as above in a form of psuedo-code, generalized for mutlipartite scenarios restricted 
to binary output and dichotomic measurements, in the manner of Eq. 



1. Accept input in the form of a linear function, for example: 

input = x(A Q ® 1) + (Ao <g> B ) + (A 1 ® B ) + (A <g> B x ) - (A x ® B x ) 



(12) 



2. Interpret statistical measures as a quantum operators, such that ® can be interpreted as the matrix Kronecker 
Product. 

Z = input : Party[fc] Apparatus[n] -> cos(9 k )a x + (-1)" sm{8 k )a y (13) 

3. Obtain the characteristic polynomial of the operator Z and all of its non-constant derivatives. 



Dq = CharacteristicPolynomial (Z) with respect to m 

(2 

= m 
3D 



(14) 



(2 K \ (2 K — 1) (2 K —2) 

TTl ' + [iirrv >+p 3 m y 1 ... + U( 2 k-i) m + P(2 K ) 



A+i = 



dm 



D 



final 



= D 



(2*--l) 



= (2 K - 1)! m 



4. Subsequent to algebraic simplification we rewrite the all expressions in terms cosine exclusively. This is prelim- 
inary to the next step where we use the cosine as a proxy for angle since we are only considering < 6 < 7t. 



Z = Z : sin 9 k -> ^1 - cos 2 6 k (15) 

5. The complete Conditions are formed by the logical union, ie AND, of all {i : .. 2 K — 1} positivity conditions. 
Conditions = Uj (V CO s0 fe , n A > 0) = V cos e fc „ Uj (A > 0) (16) 
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6. The Quantum Bound is finally obtained by minimizing m subject to all Conditions: 

QB = inf (m, such that Conditions) (17) 



Example Intermediate and Final Results from the Algorithm 

Should the reader wish to practice the algorithm as expressed above we reproduce below the positivity conditions 
and final Quantum Bound for the example input of Eq. (TTZj) . They are as follows: for all values of cos(9a) and for all 
values of cos(9b), interpreting sin^ — > yl — cos 2 9 k the characteristic polynomial variable m must satisfy 

to 4 - m 2 (2a; 2 + 8) + (x 2 - 4) 2 + 64 (cos 2 9 A ) (sin 2 9 A ) (x 2 - 4 cos 2 9 B ) (sin 2 9 B ) > 
m 3 - m (x 2 + 4) > 
3m 2 - (x 2 + 4) > 
m > 

One can readily demonstrate that the infimum of m satisfying all the above conditions is 

f|x| + 2 for|x|>2 
[V2x 2 + 8 for \x\ < 2 

as per QB 2 < - 8 ' 1 in Table I in the main text. 



APPENDIX B: Physical Symmetries 



In the bipartite and dichotomic EPR-Bell scenario [lH3j there are exactly seven independent physically motivated 
symmetries. The first symmetry is the relabeling of the parties, Alice and Bob. The next two come from the permuting 
the indices of each party's measurement apparatuses. The remaining four come from permuting the binary outputs 
any of the measurements, a redefinition of success and failure in a sense. For statistics defined the eight parameters 
of Eq. ([2]) the fundamental symmetries are given below. 

For K parties, each having access to N choices of experimental apparatus, each experiment yielding one of D 
possible outcomes, the physical symmetries are as follows: There are K\ permutations of party labeling-scheme, N\ 
permutations of the apparatus labeling scheme for each of the K parties, and D\ output labcling-schemes for each of 
the K x N apparatus settings. The total number of equivalencies, therefore, is K\ N\ K D\ KN . The total number of 
fundamental symmetries is composed of the identity, A"! — 1 party permutations, K (N\ — 1) apparatus permutations, 
and KN (D\ — 1) output permutations. 



TABLE IV. The 7 non-identity symmetries pertinent to EPR-Bell scenario. These symmetries form a group with 2 7 = 128 
equivalencies. 
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Permutations of the output 
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